Tidal effects and the Proximity decay of nuclei 
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We examine the decay of the 3.03 MeV state of 8 Be evaporated from an excited projectile-like 
fragment following a peripheral heavy-ion collision. The relative energy of the daughter a particles 
exhibits a dependence on the decay angle of the 8 Be*, indicative of a tidal effect. Comparison of the 
measured tidal effect with a purely Coulomb model suggests the influence of a measurable nuclear 
proximity interaction. 
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Aggregation of clusters in a dilute medium is a process 
that impacts a wide range of physical phenomena from 
the formation of galactic structure [? ] to formation of 
Van der Waals clusters in a low density gas[? ]. This 
aggregation can involve the delicate interplay of elemen- 
tary forces that results in a frustrated system. One such 
example is the formation of pasta nuclei in the crust of a 
neutron star0,0. For smaller systems, the phenomenon 
of alpha clustering is important both in low density nu- 
clear matter [jjj], as well as in light nucleiQ. In the case 
of heavy nuclei, cluster aggregation is manifested in the 
spontaneous phenomenon of clusters, from the common 
process of a decay to the emission of more exotic clusters 
such as 14 C The reduction of density near the nuclear 
surface, allows formation of a particles, or other clusters, 
and their emission from either ground-state or modestly 
excited nuclei. Cluster emission, thus primarily probes 
the surface properties of the emitting nucleus 6]. Our 
present understanding of cluster emission is largely based 
upon the yields, kinetic energy spectra, and angular dis- 
tributions of emitted clusters - all of which are well de- 
scribed within a statistical transition-state formalism^. 
In this Letter we present for the first time evidence for the 
modification of cluster emission by interaction with the 
nuclear surface. We probe the interaction of the nuclear 
surface with the emitted cluster by using resonance spec- 
troscopy to examine the emission of 8 Be* and specifically 
explore how its decay is impacted by the tidal effect @. 

Charged-particles produced in the reaction 
114 Cd+ 92 Mo at E/A=50 MeV were detected in an 
exclusive Air setup. To focus on evaporated fragments, 
we selected peripheral collisions through detection of 
forward-moving projectile-like fragments (PLFs) with 
10<Z<48, in the angular range 2.r<9 tab <4.2° with 
A(9 iab «O.13°0. This PLF is the decay residue of 
the excited primary projectile- like fragment (PLF*) 
formed by the collision. Light-charged-particles and 



fragments with Z<9 were isotopically identified 9] in the 
angular range 7°<9 lab <58° with the silicon-strip array 
LASSA [H El. Each of the 

nine telescopes in this 
array consisted of a stack of three elements, two 5cm x 
5cm silicon strip detectors (Si(IP)) backed by a 2 x 2 
arrangement of CsI(Tl) crystals each with photo-diode 
readout. Each telescope was segmented into 16x16 
orthogonal strips, resulting in good angular resolution 
(A6> iab «0.43°). The 3 x 3 arrangement of the LASSA 
telescopes was centered at a polar angle 9 lab =32° with 
respect to the beam axis. The identification threshold 
of LASSA is 2 MeV/ A for a particles. Light-charged 
particles emitted at other angles were detected by the 
Miniball/Miniwall Air phoswich array. For the analysis 
that follows events were selected with 15<Zp£p<46, 
Vp_L_F>8.0cm/ns and the multiplicity of particles in 
LASSA, Ni J 455A=2,3. Using the measured emitted 
particles and the assumption of isotropy, the Z, A, and 
velocity of the PLF* was calculated \L2\. Associated with 
these events is a PLF* with a most probable atomic 
number of «30. 

The kinetic energy spectra of a particles and Be nuclei 
emitted in the angular range #plp»<80° are presented 
in Fig. ^ These spectra are reasonably well described 
by a Maxwcll-Boltzmann distribution supporting the pic- 
ture of a largely evaporative process [l2j. In comparison 
to the a particle spectra, the spectra of Be fragments 
exhibit more pronounced tails indicative of a higher ini- 
tial temperature 01 ■ The average kinetic energy of 7 Be 
is larger than that of 9,10 Be in agreement with previ- 
ous work suggestive of sequential decay of excited pri- 
mary fragments as they propagate away from the emit- 
ting nucleus 9] . For the events presented, damping in the 
collision process produces a PLF* with a most probable 
excitation energy of 3 MeV/ A and a maximum excita- 
tion of 4 MeV/A. The deduced excitation is consistent 
with the "temperature" associated with the a particle 
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FIG. 1: Kinetic energy spectra of a particles and 7,9,10 Be 
fragments in the PLF* frame observed with 6plf* <80°. 
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FIG. 2: Panel a: Relative kinetic energy distributions for 
a- a pairs. Panel b: Comparison of the results of a non- 
resonant final-state interaction Monte Carlo calculation with 
the mixed-event background. 



kinetic energy spectrajT^]- The de-excitation cascade of 
the PLF* involves the emission of not only nucleons and 
ground state fragments, but also excited fragments, some 
of which are unstable against particle emission. 

The decay of short-lived resonant states such as 8 Be* 
can be examined by constructing the relative energy spec- 
trum of the daughter products. Shown in Fig. |^ is 
the relative energy spectrum between two a particles 
for events in which either two or three a particles were 
detected in LASSA with 9 atPLF »<100°. Solid symbols 



depict the experimental data which exhibits a peak at 
E re ;«3 MeV. The overall shape of the E re ; spectrum is 
affected at low E re ; (< 2 MeV) by the finite angular ac- 
ceptance of the LASSA CsI(Tl) detectors, as two parti- 
cles entering the same CsI(Tl) crystal are not resolved. 
The 8 Be ground-state is consequently not observed in 
this kinematic regime. For larger E re ; (>5 MeV), the de- 
crease in yield is impacted by the geometric acceptance 
of LASSA. The observed relative energy distribution can 
be understood as having two primary components: res- 
onant decay of 8 Be* and non-resonant a emission. In 
order to assess the non-resonant contribution, which acts 
as a "background" for the resonant decay, we performed a 
mixed-event analysis. Two a particles were selected from 
different events and their relative energy was calculated. 
The resulting relative energy spectrum was normalized 
in the interval 14<E re ;<50 MeV (where the resonant de- 
cay contribution is expected to be insignificant) and is 
displayed as open squares in Fig. 12k 1 -- While the mixed- 
event background has the general shape of the observed 
relative energy spectrum, there is an excess yield in the 
observed data centered at «3 MeV. This excess originates 
from the decay of the first excited state of 8 Be at 3.03 
MeV with an intrinsic width of 1.5 MeV. Over-prediction 
of the yield for E re ;<1.5 MeV is understandable as the 
mixed-event pairs do not experience Coulomb repulsion. 
One of the drawbacks of the mixed-event analysis is that 
the events sampled span a broad distribution of PLF* 
velocities. To select the conditions that best describe 
the non-resonant a emission, we restricted the differ- 
ence in PLF* velocity, AVp^p*, for the events in the 
mixed-event background. The mixed-event background 
corresponding to AV plp*<0.6 cm/ns is shown as open 
triangles in Fig. 01. This background provides a better 
description of the data for E re ;<1.5 MeV while it ex- 
hibits a larger yield at E re ;«3 MeV. To determine which 
value of AVplf* provided the best reproduction of the 
non-resonant decay, we calculated the \ 2 between the 
mixed events and the experimental data in the interval 
6<E re ;<50 MeV. The minimum x 2 /^=l-l-l corresponds 
to AVp LF *= 0.4-0.6cm/ns. 

In order to better understand the "background" , we 
have modeled the non-resonant final-state interaction 
with a Monte Carlo Coulomb trajectory calculation. In 
this model, which we refer to as MC-FSI, two a particles 
are isotropically emitted in sequential fashion from the 
surface of an excited nucleus following Lambert emission, 
while accounting for all recoil and Coulomb interactions. 
The Zplf*, A-plf*, Vplf*, and Oplf* are taken from 
the experimental data while the time between emissions 
is taken to be exponential with a mean time r. Fol- 
lowing Coulomb propagation all particles were filtered 
for the experimental acceptance, angular resolution, and 
thresholds. Calculations with r=100, 300, and 500 fm/c 
are shown by the lines in Fig. along with the ex- 
perimental data and mixed-event background. With de- 
creasing mean emission time, suppression of yield for low 
values of E re ; is observed. While the calculation with 
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FIG. 3: Difference relative energy spectra of observed a pairs 
for different decay angles, /3, relative to the mixed-event back- 
ground. 



r=100 fm/c is clearly inconsistent with the experimen- 
tal data, r>300 fm/c provides a reasonable description 
of the mixed-event background for E re ;>2 MeV. A de- 
duced mean emission time of this magnitude associated 
with an excitation energy of 3 MeV/A is consistent with 
previous work^^. For smaller values of E re ;, mutual 
Coulomb repulsion of the two a particles results in a 
suppression of yield at small E re ; as compared to the 
mixed-event data. The comparison of the mixed-event 
background with the Monte Carlo final state calculation 
indicates that for all but the smallest values of E re ; the 
mixed-event background provides a reasonable descrip- 
tion of the non-resonant contribution. 

The first excited state of 8 Be has an intrinsic width of 
1.5 MeVfLU corresponding to a mean lifetime of « 131 
fm/c. For such a short lifetime the likelihood that this 
state decays in the vicinity of the emitting nucleus is 
significant. Consequently, as the 8 Be* decays into two 
a particles, its increasing quadrupole moment interacts 
with the gradient of the Coulomb field. For a pairs which 
decay orthogonal to the emission direction, the field gra- 
dient provided by the emitting nucleus acts to increase 
their relative energy. In contrast, a pairs which decay 
along the emission direction experience a reduced rela- 
tive energy due to the larger acceleration of the nearer 
a particle Q. In order to examine this Coulomb tidal 
effect, we constructed the difference relative energy spec- 
tra of the observed a pairs relative to the mixed-event 
background. These spectra were constructed for differ- 
ent decay angles, /?, calculated as the angle between the 
relative velocity of the a-a pair and the center-of-mass 
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FIG. 4: Dependence of <E re ;>, associated with the two a 
particles arising from the decay of first excited state of 8 Be, 
on the average decay angle, </?>. Lines represent the results 
of Monte Carlo calculations for resonance decay. 



momentum of the pair. Depicted as the solid symbols 
in Fig. |21 is the case using the mixed-event background 
with AVplf*<0.6 cm/ns. For reference, the difference 
spectra with no restriction on AVpLp, are also shown 
as open symbols. Clearly the AVplf* restriction docs 
not introduce a major change in the difference relative 
energy spectra. The prominent feature of the difference 
spectra is the peak at «3 MeV which is evident for all 
p. The negative dip for E re ;<«1.5 MeV in the difference 
spectra is due to the over-prediction of the yield by the 
mixed-event background. Moreover, it is clearly evident 
that the peak in the difference spectrum for 0</3<45° 
is shifted to lower values of E r ei as compared to larger 
values of [3. To quantitatively extract the dependence of 
<E re i> on for the 3.03 MeV state, we integrated the 
difference spectra over the region indicated by the verti- 
cal lines. These integration limits were selected by con- 
sidering the intrinsic width of the 3.03 MeV state and the 
shape of the measured spectra. The extracted <E re ;> is 
largely insensitive to any reasonable choice of integration 
limits. Using these difference spectra and accounting for 
efficiency, we also extracted the total yield for the 3.03 
MeV state. Based upon the yield of measured 7 ' 9 Be iso- 
topes we estimate the yield for the 3.03 MeV state. The 
extracted yield agrees with this expectation to within 
15%. 

The dependence of the average relative energy, 
<E re j>, on the decay angle j3 is shown in Fig. 0] as solid 
symbols. A clear manifestation of the tidal effect is ob- 
served as <E re ;> increases with increasing (3. The mag- 
nitude of the observed change in <E re ;> is «0.7 MeV, a 
relative change of «25%. For reference, the arrow indi- 
cates the decay energy corresponding to the intrinsic en- 
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ergy of the state and the Q- value of the decay. The error 
bars shown, while dominated by the measurement statis- 
tics, include the uncertainty associated with the integra- 
tion limits. To demonstrate the sensitivity of <E re ;> to 
different mixed-event backgrounds we also show as open 
symbols the dependence of <E re ;> on (3 for no restric- 
tion on AVplf*- The same overall trend of <E re ;> on 
(3 is observed for the different mixed-event conditions. 

To understand this observed trend quantitatively, we 
simulated the decay of a 8 Be* in a simple Monte Carlo 
model, MC-Res. In this model the 8 Be* is emitted 
isotropically from the surface of the PLF*. The prop- 
erties of the PLF* were sampled in the same manner as 
in the MC-FSI calculations. The lifetime of the emitted 
8 Be* was chosen consistent with an exponential proba- 
bility distribution P(t)=exp(-t/r), reflecting first order 
kinetics. The mean lifetime, r, was taken to be 131 fm/c 
determined by the intrinsic state width. The initial ki- 
netic energy of the 8 Be* is taken to be exponential with 
a slope parameter of 7.5 MeV, consistent with the exper- 
imental data shown in Fig. ^ The 8 Be* propagates in 
the Coulomb field of the emitting nucleus until it decays. 
At the moment of decay the 8 Be* is replaced with two 
a particles with an inter-a separation distance (scission 
configuration) given by R a ~ a = 5.81 fm in accordance 
with svstematics|l5|. Selecting a smaller quadrupole mo- 
ment, namely Rq,_ q = 4 fm, makes a negligible difference 
in the final results. The decay angle, (3, of the two a 
system with respect to the emission direction is taken 
to be isotropic namely the effect of the Coulomb field 
in orienting the decaying 8 Be* is neglected. Following 
decay, the two a particles are propagated along trajecto- 
ries corresponding to both the Coulomb repulsion from 
the PLF* and their mutual Coulomb repulsion. Energy 
and momentum are conserved at all stages of the calcu- 
lation. Particles are subsequently filtered by the detector 
acceptance, angular resolution, and thresholds. 

The results of this schematic calculation are shown 
as the dashed line in Fig. 21 The angular dependence 
of <E re i>, namely the magnitude of the tidal effect, is 
clearly overestimated by the purely Coulomb calculation. 
One possible reason for this difference is the interaction 
of the emitted 8 Be* with the nuclear surface. This in- 
teraction is particularly important as the most probable 
initial velocity of the 8 Be* is close to zero, and conse- 
quently the 8 Be* spends a significant time within the 
range of the proximity interaction. This is in marked 
contrast to the tidal decay following projectile breakup 
reactions|l6j. The attractive surface interaction can im- 



pact the decay of the evaporated 8 Be* in two ways. It 
can stabilizes the 8 Be* while it is in the vicinity of the 
emitting nucleus. In addition, the attractive potential 
may produce a nuclear tidal effect with the opposite sign- 
as compared to the Coulomb tidal effect. Both the nu- 
clear tidal effect and the increased stability of the 8 Be* 
will result in a decreased dependence of <E re ;> on (3. 
To assess whether nuclear attraction could be respon- 
sible for the difference between the observed tidal effect 
and the Coulomb calculation, we have introduced a delay 
into the decay time distribution such that P(t)=exp(-(t- 
^•deiay) / t) for t>tdeiay This artificial delay mimics the 
increased stability of the 8 Be* due to the nuclear interac- 
tion. With increasing tdeiay, the dependence of <E re ;> 
on (3 decreases as evident in Fig.0] For td e / ay =200 fm/c, 
represented by the solid line, one observes reasonable 
agreement with the observed tidal effect. Without any 
delay the lifetime of 131 fm/c corresponds to a decay dis- 
tance of «17 fm. For t ( ; e ; a j / =200 fm/c, the 8 Be* travels 
on average an additional 4 fm before decaying. From this 
we conclude that stabilization of the 8 Be* by the emitting 
nucleus for a distance of «4 fm is sufficient to reproduce 
the magnitude of the observed tidal effect. A more realis- 
tic model of the nuclear interaction is necessary in order 
to extract more quantitative results. 

In summary, we have measured the dependence of 
<E re j> on decay angle for the first excited state of 8 Bc. 
This trend can be understood as a tidal effect in which 
the decaying cluster interacts with the gradient of an 
external field provided by the emitting nucleus. Calcu- 
lations with a Coulomb model over-predict the observed 
tidal effect. This over-prediction suggests that the emit- 
ting nucleus influences the cluster decay through the nu- 
clear interaction. The attractive nuclear potential of the 
emitting nucleus acts to both stabilize the excited clus- 
ter and induces a nuclear tidal effect which opposes the 
Coulomb tidal effect. 
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